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O ■ Abstract 

We calculate the time dependence of B° — > / decays (and analogously D — > / decays) 
and use this result to obtain expressions for the time dependence of: an untagged sample 
(B° — > / and B — > / decays combined), a sample containing charge-conjugate final states 
(B° — ► / and B° — > / decays combined), and a sample containing all four decay modes 
\q ' together. For simplicity we assume CP violating effects are negligible. 
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The time dependence of B° — » / decays, in which / is not a CP eigenstate, is 
not purely exponential due to the presence of B°-B° mixing. This mixing arises 
due to either a mass difference Am or a decay-width difference A7 between the 
mass eigenstates of the B°-B system. This note calculates an expression for the 
time dependence in the presence of such mixing. It is assumed that CP violation, 
^h. if present, occurs at a negligible level. This assumption is well-motivated in the 

D° system, for which our results can be applied. In particular, our results can be 
applied to the decay D° — > K~ir + ||]. 

There are two strong eigenstates, B° and B° , which are not eigenstates of the 
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5-H 



full Hamiltonian due to the weak interaction. Thus, they most generally evolve 



according to the Schrodinger equation: 
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where the 2x2 Hermitian matrices M and T represent transition amplitudes to 
virtual (off-mass-shell) and real (on-mass-shell) states, respectively. Diagonalizing 
M — (i/2)r yields the physical states 

\B S ) = p\B ) + q \B°) 

\B L ) = p\B°)- q \B°) (2) 

where 



q _ l M* 2 -(i/2)T* 12 ^ Mh _ ™ 



p V M 12 -(i/2)r 12 V M : 



e lI *M . (3) 
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In Eq. (0) we've used the notation B L /B S in analogy with the K L /K S system, 1 
while in Eq. (||) we've made the good approximation T 12 <C M 12 f|]. The time 
evolution of the physical states is: 

\B s (t)) = \B s )e- irs/2+im s^ 

\B L (t)) = \B L ) e -^L/2+im L )t _ ( 4 ) 



Inverting Eq. (0) gives: 



\B°) = ^(\B S ) + \B L )) 

\B°) = Y q (\B s )-\B L )) , (5) 



and inserting the time dependences (||) into Eqs. (||) gives: 
\ B °(t)) = l{| J B 5 } e -( r s/2+- s )t + | jBi ) e -(r i /2+ imi ) i | 

= e -(r/2+im)t | cogh [(A 7 /4 + iAm/2)t] \B°) + (^j sinh [(A 7 /4 + iAm/2)t] 
\B°(t)) = l{|5 5 )e-( r s/ 2 +™s)* -\B L )e^l 2+im ^) 

= e -(r/2+im)< j^j s inh[(A 7 /4 + iAm/2)t] \B°) + cosh [(A 7 /4 + iAm/2)t] \B°)} 

where m = (m L + m s )/2, T = (T L + T 5 )/2, Am = m L — m s , and A 7 = T L — T s . 
In the limit A 7 — > 0, one recovers the expressions which follow Eq. (5) in Ref. g. 
The states above lead to decay amplitudes: 

(f\H\B°(t)) = e ~^ /2+m) 1 jcosh [(A 7 /4 + iAm/2)t] A f + (^j sinh [(A 7 /4 + iAm/2)t] A f 
(f\H\B°(t)} = e - {T/2+m] 1 | sinh [(A 7 /4 + »Am/2)i] A f + cosh [(A 7 /4 + iAm/2)t] A f 



where we've defined the amplitudes for the pure B° and B° states as: 

A f = (f\H\B°) Af = (f\H\B°) . (6) 

Squaring the decay amplitudes gives the decay rates: 

\(f\H\B°(t))\ 2 = \A f \ 2 e~ Tt {|cosh(A 7 /4 + iAm/2)t | 2 + 

|A| 2 |sinh(A 7 /4 + iAm/2)t | 2 + 

(A*)cosh(A 7 /4 + iAm/2)i [sinh(A 7 /4 + iAm/2)t ]* + c.c.j 

(7) 

Although \t Bl - r B J < \t k ^ -T Kg \. 
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\{f\H\B°{t))\ 2 = \A f \ 2 e- rt {|cosh(A 7 /4 + iAm/2)t | 2 + 

|A| 2 |sinh(A7/4 + iAm/2)t | 2 + 

(A*) cosh(A7/4 + iAm/2)t [sinh(A7/4 + iAm/2)t }* + c.cj , 

(8) 

where we've defined the parameters A = (q/p)(Aj/Aj) and A = (p/q)(Af/ Af). To 
evaluate these expressions, note that: 



cosh(A7/4 + iAm/2)t\ 



| cosh(A7/4)t cos(Am/2)t + i sinh(A7/4)t sin(Am/2)t[ 2 
cosh 2 (A7/4)t cos 2 (Am/2)t + sinh 2 (A7/4)t sin 2 (Am/2)t 

cosh 2 (A7/4)t - sin 2 (Am/2)i 
1 
2 



cosh(A7/2)t + cos(Am)t 



(9) 



Similarly, 

| sinh(A7/4 + iAm/2)t| 2 



| sinh(A7/4)t cos(Am/2)t + i cosh(A7/4)t sin(Am/2)t[ 2 
sinh 2 (A7/4)t cos 2 (Am/2)i + cosh 2 (A7/4)t sin 2 (Am/2)i 
cosh 2 (A7/4)t - cos 2 (Am/2)i 
V 



Finally, 



cosh(A7/2)t — cos(Am)i 



cosh(A7/4 + iAm/2)t [sinh(A7/4 + iAm/2)t]* 



(10) 



[cosh(A7/4)i cos(Am/2)t + i sinh(A7/4)t sin(Am/2)t] x 

[sinh(A7/4)i cos(Am/2)t - i cosh(A7/4)t sin(Am/2)i] 



sinh(A7/2)t cos 2 (Am/2)t + sinh(A7/2)t sin 2 (Am/2)t + 

isinh 2 (A7/4)tsin(Am)i — i cosh 2 (A7/4)t sin(Am)i 
sinh(A7/2)t — isin(Am)i 



(11) 



Inserting Eqs. @, @, and flllj) into Eqs. (0) and (§) gives: 

\(f\H\B°(t))\ 2 = ^^^e- r *{[cosh(A7/2)t + cos(Am)t] + 

|A| 2 [cosh(A7/2)t - cos(Am)t] + 
(A*) [sinh(A7/2)t -ism(Am)t } 
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(A) [sinh(A7/2)t + isin(Am)t ]} 

= (^j-^j e~ Tt {(l + |A| 2 ) cosh(A 7 /2)t + (1 - |A| 2 ) cos(Am)t + 

(A + A*) sinh(A 7 /2)t + i(X - A*) sm(Am)t } 

(12) 

\(f\H\B°{t)}\ 2 = (^j-^j e- r *{[cosh(A 7 /2)t + cos(Am)t] + 

|A| 2 [cosh(A7/2)t - cos(Am)t] + 

(A*) [sinh(A7/2)t -zsin(Am)t] + 

(A) [sinh(A7/2)t + i sm(Am)t ]} 
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-Tt 



{(1 + |A| 2 ) cosh(A 7 /2)t + (1 - |A| 2 ) cos(Am)t + 

(A + A*) sinh(A7/2)t + i(X - A*) sin(Am)t } . 

(13) 

If the lifetime distribution for a final state / is constructed from a sample of 
untagged decays, and there are equal numbers of B° and B° mesons produced, 
then the time dependence of the decays will be the sum of Eqs. (0) and (U). If 
\Af \ <C \AA (e.g., if B° — > f is Cabibbo-favored and B° — > / is doubly-Cabibbo- 
suppressed), then one would measure: 

- Q)e--{2|^| 2 cosh(A7/2)t + 

A}A f [(q/p) + (p/q)*} sinh(A 7 /2)t + A f A} [(q/p)* + (p/q)] sinh(A 7 /2)t + 
A* f A f [(q/p) - (p/q)*]ism(Am)t + A f A* f [(p/q) - (q/p)*]ism(Am)t } 

^ e~ Tt {2|^l / | 2 cosh(A7/2)t + \A}A f (2q/p) + A f A* f (2q/p)*^ sinh(A 7 /2)t } 

|^ / | 2 e~ r '{cosh(A 7 /2)t + (A + A*) sinh(A 7 /2)t } , (14) 



where we've used the fact that \q/p\ 2 = 1. This result [and Eq. ( |l5| ) below] was 
previously obtained by Dunietz Hj. 

If the lifetime distribution includes final states / and / combined together, then 
we must also consider the decay rates \(f\H\B° (t))\ 2 and \(f\H\B°(t)}\ 2 . These rates 



are equivalent to Eqs. (|12|) and (|13|) , respectively, with the interchange q/p <-> p/q 
(since \Aj\ = \Aj\ and \Af \ = \Aj\ by CP conservation). The analog of Eq. ( |l4| ) 
for the final state / is then: 
dN _ _ 

H/ « |^ / | 2 e- r *{cosh(A 7 /2)t + ( K + K *)sinh(A 7 /2)i}, (15) 
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where « = (p/q)(A f -/A f -) = (q/p)* (A f /A f )e ie . Thus, 
dN _o 

(g ~ 2|^ / | 2 e~ r '{cosh(A 7 /2)t + i?e(A + k) sinh(A 7 /2)t } . 

(16) 

For this expression to be valid the / and / final states must have equal acceptances, 
or else the decays must be corrected for acceptance before the lifetime distributions 
are combined. 

If the B° and B mesons are produced in a fixed-target experiment, then their 
yields will usually be different due to there being unequal numbers of d and d quarks 
in the initial state. In this case Eq. ( |I4| ) does not apply, and it is more convenient 
to consider B° and B° samples separately. Such experiments typically combine 
together the final states / and / (when not studying CP asymmetries), and in this 
case: 

dA W(/+/) _ . |2 ft 



dt 



dt 



A f \ 2 e- rt {(l + |A| 2 ) cosh(A 7 /2)t + (17) 
Re(X + «*) sinh(A 7 /2)t - Im(X + k*) sin(Am)t| 

\A f \ 2 e- ft { (l + |A| 2 ) cosh(A 7 /2)t + (18) 
Re(X* + k) sinh(A 7 /2)t - Im(X* + k) sin(Am)t| . 



Adding together Eqs. (|T^) and ((T^) (corresponding to an untagged sample having 
equal numbers of B° and B ) recovers Eq. (|l6[) . 
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